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Abstract. - We analyze the properties of the energy landscape of finite-size fully connected 
spin-glass models with a discontinuous transition. In the thermodynamic limit the equilibrium 
properties in the high temperature phase are described by the schematic Mode Coupling Theory 
of super-cooled liquids. We show that finite-size fully connected spin-glass models do exhibit 
properties typical of Lennard- Jones systems when both are near the critical glass transition, 
where thermodynamics is ruled by energy minima distribution. Our study opens the way to 
consider activated processes in real glasses through finite-size corrections (i.e. calculations 
beyond the saddle point approximation) in mean-field spin-glass models. 



In recent years a significant effort has been devoted to the understanding of glass-forming 
systems. Recent theoretical and numerical results clearly show that the slowing down of 
the dynamics is strongly connected to the potential energy landscape. The trajectory of the 
representative point in the configuration space can be viewed as a path in a multidimensional 
potential energy surface. The dynamics is therefore strongly influenced by the topography of 
the potential energy landscape: local minima, barriers heights, basin of attraction an other 
topological properties all influence the dynamics. 

The potential energy surface of a super-cooled liquid contains a large number of local 
minima, called inherent structures (IS) [Q, each surrounded by a basin defined as the set of all 
configurations that a local energy minimization maps onto the IS contained within it. In this 
picture the time evolution of the system can be divided into intra-basin and inter-basin motion. 
Transition from one basin to another are expected to occur differently as the temperature is 
varied. In particular, when the temperature is lowered down to the order of the critical Mode 
Coupling Theory (MCT) temperature Tmct the two motions become well separated in time 
and the relaxation dynamics is dominated by slow thermally activated crossing of potential 
energy barriers between different basins ||]. 
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The essential features of MCT for glass-forming systems are also common to some fully 
connected spin glass models, called mean-field p-spin glass with p > 2 ||^|]. In the thermody- 
namic limit the high-temperature phase, (paramagnetic in the spin-glass language and liquid 
in glass language) , is described by the schematic mode MCT for super-cooled liquids . As 
a consequence at the critical temperature, called Tjj in p-spin language and Tmct in MCT 
language, an ergodic to non-ergodic transition takes place. In mean-field models the barriers 
separating different ergodic componets are proportional to the system size, hence infinite in 
the thermodynamic limit, and at Tjj the relaxation time diverges. We note that despite this 
strong ergodicity beaking the dynamical behaviour is far from being trivial [Q ^| . In real sys- 
tems the barriers are of finite height and the transition to a glassy state appears at T g < Tjj, 
the glass transition temperature, where the typical activation time over barriers is of the same 
order as the observation time. 

Despite these differences mean-field models, having the clear advantage of being analyti- 
cally tractable, have been used to study the phase space structure of glassy systems, especially 
between the dynamical temperature Td and the static temperature T c (Kauzmann tempera- 
ture Tk in glass language) where a real thermodynamic phase transition driven by the collapse 
of the configurational entropy takes place. The main drawback is that, since activated process 
cannot be captured by mean-field models, the picture that emerges is not complete. Let us 
remark also that, despite the large amount of analytical work devoted to the study of the 
static as well as dynamical properties in the N — > oo limit, much less is known concerning the 
finite N behavior. 

In this Letter we investigate numerically finite-size fully connected spin-glass models, where 
activated process are present. We find strong indications that, once activated process are 
allowed, mean-field spin glass models with a discontinuous transition |1 do exhibit properties 
similar to the ones recently observed in simulations of Lennard- Jones systems near Tmct 
|t]-|1(|, making these systems highly valuable for studying the glass transition not only at 
MCT level but also in the activated regime. 

Here we shall focus on the equilibrium properties of the Ising-spin Random Orthogonal 
Model (ROM) JTT|,|l^ ] which shows a discontinuous spin-glass transition identical to that found 
in p-spin models ]13|. Similar results are, however, obtained using other p-spin-like models, as 
for example the Bernasconi Model and the Ising p-spin model. Consequently, our conclusions 
and results should be generally valid for mean-field spin glasses with a discontinuous transition. 
The advantage of the ROM lies in its interaction term which is two-body, at difference with 
the p-body interaction of p-spin models, making the simulations much faster. Moreover, the 
ROM has a very strong freezing glass transition Jll]] making it an extremely good microscopic 
realization of the Random Energy Model [Q . 

The model is defined by the Hamiltonian 

H=-2'£j ij v i * j (1) 

ij 

where <7j = ±1 are N Ising spin variables, and Jy is a TV x N random symmetric orthogonal 
matrix with Ju — 0. Numerical simulations are performed using the Monte Carlo method 
with the Glauber algorithm. For N — > oo this model has the same thermodynamic properties 
of the p-spin model: a dynamical transition at Td = 0.536, with threshold energy per spin 
eth = Eth/N — —1.87, and a static transition at T c = 0.25, with critical energy per spin 
e lrsb = -1.936 |llj,p. 

The two transitions can be understood from a "geometrical" point of view. The TAP 



analysis of these models [12,15 reveals that the phase space visited is composed by an expo- 



nentially large (in N) number of different basins, separated by infinitely large (for N — > oo) 
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barriers. Different basins are unambiguously labelled by the value of the energy density e at 
T = 0. Each TAP solution describes the thermodynamics within the basins labelled by e and 
at T = coincide with the potential energy local minima, i.e. the IS of the system. The 
dynamical transition is associated with IS with the largest basin of attraction for A — > oo, 
while the static transition with IS with the lowest accessible free energy 

(H Ell- 



in the mean-field limit, the allowed values of e are between e\ rs b and eth- Solutions with 
e larger than eth are unstable Jl7] ], while solutions with e smaller than ei rs b have negligible 
statistical weight Jl2],|l5| ■ The IS with energy equal to the threshold energy eth can be identified 
with the less stable solutions. These IS have the largest (exponentially with N) statistical 
weigth so that a rapid cooling from an equilibrium state at T > Try down to T = will, with 
probability one for A" — > oo, drives the system to IS with e = e t h- For finite A~ the scenario is 
different since not only IS with e < eth acquire statistical weigth, but solutions with e > eth 



and few negative directions 17 may become stable, simply because there are not enough 
degrees of freedom to hit them. 

To analyze the thermodynamics of finite systems we follow the ideas of of Stillinger and 



Weber 18 and decompose the partition sum into a sum over different IS and a sum within 
each basin. Collecting all IS with the same energy e, denoting with exp[Ns c (E)] de the number 
of IS with energy between e and e + de, and shifting the energy of each basin with that of the 
associated IS, the partition sum can be written as jl8| 

Z N (T)~ J rfeexp A [-(3e+ s c (e) - Pf(J3,e)] (2) 

where f(f3, e) can be seen as the free energy density of the system when confined in one of 
the basin associated with IS of energy e. The function s c (e) is the configurational entropy 
density also called complexity (the same derivation can be also obtained within the replica 
approach Jl9| ) . From eq. (^) we easily obtain the probability that an equilibrium configuration 
at temperature T = 1/(3 lies in a basin associated with IS of energy between e and e + de: 

P N (e, T) = exp AT \-(3e + s c (e) - (3f(P, e)] /Z N [T). (3) 

Taking the N oo limit we recover the mean-field results Here we keep A finite. 

From the partition function we can compute the average internal energy density. This is 
given by u(T) = (e + d((3f)/d[3) — (e(T)) + (Ae(T)), where the first term is the average energy 
of the IS relevant for the thermodynamics and the second is the contribution coming out from 
fluctuations within the IS. The average is taken with the weight In the limit A" — > oo the 
only relevant IS are those with e = eth, and limAr^oo (e(T)) = eth for any T > Td- To measure 
(e(T)} for finite A^ we consider the following experiment. First we equilibrate the system at 
a given temperature T , then we instantaneously quench it at T = 0. This is obtained by 
decreasing the energy along the steepest descent path. In this way we can identify the energy 
of the IS. The experiment is repeated several times starting from uncorrelatcd equilibrium 
configurations at T and the average IS energy is computed. In figure [l] we report (e(T)} as a 
function of temperature T for system sizes A^ = 48, 300 and 1000. As expected, as A* increases 
(e(T)} tends to eth- From the numerical data we found that the plateau energy, approaches 
e t h with the power law (e p / Qteau ) — e t h ~ A~ a2 , see inset Fig. [j] Note that since A^ is finite 
we can equilibrate the system also below Try, down to the glassy transition T g (N), ~ 0.35 for 
A^ = 48 and ~ 0.5 for N — 300, below which the system falls out of equilibrium p0| . 

The figure shows that for finite A^ and T not too close to Td the thermodynamics is 
dominated by IS with e > eth- This is clearly seen from the (equilibrium) probability dis- 
tribution of e since it is centered about (e(T)) indicating that IS with e ~ (e(T)) have the 
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Fig. 1 - Temperature dependence of (e(T)> for N = 48 (square), N = 300 (circle) and N = 1000 
(triangle). The average is over 10 3 different equilibrium configurations at temperature T. The hori- 
zontal line is the N — + oo limit. The arrows indicate the critical temperatures To and T c (see text). 
The dotted line is the curve obtained from the configurational entropy for large TV. In the inset it 
is shown the size dependence of (e(T = 3)) as a function N~ a with a = 0.2 extrapolated down to 
-1.87 (triangle). 



largest statistical weight. This scenario has been also observed in numerical simulation of, 
e.g., Lennard- Jones systems 00,53,5^]. 

In the temperature range where eq. is valid, the curves lnP/v(e,T) + /3e are equal, 
except for a temperature dependent factor lnZjv(T), to s c (e) — /3/(/3,e). If we can neglect 
the energy dependence of f((3,e), then it is possible to superimpose the curves for different 
temperature. The resulting curve is, except for an unknown constant, the complexity s c (e). 
The curves obtained for system sizes N = 48 and 300 and various temperatures between 
T = 0.4 and T = 1.0 are shown in figure || (a). The data collapse is rather good for e < —1.8. 
Above the curves cannot be superimposed indicating that the energy dependence of /(/3, e) 
cannot be neglected anymore. In liquid language this is the anharmonic threshold f|,[l(]]. To 
compare the result with the analytical predictions each curve has been translated to maximize 
the overlap with the theoretical prediction for s c (e) The line is the quadratic best fit 

from which we can estimate the critical energy e c ~ —1.944, where s c (e c ) vanishes, in good 
agreement with the theoretical result ei rs b = —1.936 |0|. 

Direct consequence of f(/3,e) ~ f((3) for e < —1.8 is that in this range the partition 
function can be written as the product of an intra-basin |p3[ contribution [exp(— N(3f)] and of 
a configurational contribution which depends only on the distribution of IS energy densities. 
The system can then be considered as composed by two independent subsystems: the intra- 
basin subsystem describing the equilibrium when confined within basins, and the IS subsystem 
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Fig. 2 - (a) Configurational entropy as a function of energy. The data are from system sizes N = 48 
and N = 300, and temperatures T = 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. For each curve the unknown 
constant has been fixed to maximize the overlap between the data and the theoretical result jl^] . The 
line is the quadratic best-fit. (b) Configurational entropy density as a function of temperature. The 
line is the result from the best-fit of s c (e) while the symbols are the results from the temperature 
integration of eq. (|) for N = 48, N = 300 and N = 1000. 



describing equilibrium via activated processes between different basins. As the temperature is 
lowered and/or N increased the two processes get more separated in time and the separation 
into two subsystems becomes more and more accurate. A similar scenario is also oserved 
in numerical simulations of super-cooled liquids, such as Lennard-Jones systems, near the 
MCT transition j2|,|24|]. The form of f(/3) for the specific system can be computed studying 
the motion near the IS, for example using an harmonic approximation (^,0. However, this 
usually gives only a small corrections to thermodynamic quantities for T close to Trj P, [To| 
and we do not consider it in this Letter. 

Another important consequence of the separation into two subsystems is that in eq. ([|) 
f(J3, e) can be neglected since it cancels with the equal term coming from the denominator. 
Therefore from the knowledge of s c (e) we can easily compute the average IS energy density 
(e(T)). For large N this is given by the saddle-point estimate, see eq. (||): 

max [— /3e + s c (e)] . (4) 

e 

The result obtained using for s c (e) the quadratic best-fit, is the dotted line shown in Fig. |l|. 
The agreement with the direct numerical data is good already for A" = 300. From the form of 
(e(T)} for large N we can identify the static critical temperature T c ~ 0.3 as the temperature 
below which (e(T)) remains constant, not far from the theoretical results T c = 0.25 |l2| ]. 
To have a check of our results we have computed the configurational entropy density using a 
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Fig. 3 - Temperature dependence of (e(T)) for the ±J SK model and N = 192 (square), N = 320 
(circle) and N = 640 (triangle). The average is over 10 3 different equilibrium configuration at 
temperature T. Inset: Configuration entropy as a function of energy for the ± J SK model. The data 
are from system sizes N = 64 (circle) and N = 192 (square), and temperatures T = 0.8, 0.9, 1.0, 1.1, 
1.2, 1.3. For each curve the unknown constant has been fixed to maximize the overlap between the 
data and the theoretical result (continuous line), pa]. 



different approach |9|] . When the IS subsystem is in thermal equilibrium then the temperature 
dependence of the configurational entropy density can be evaluated from the thermodynamic 
relation 

d ^(T) 1 

d(e(T)) T 1 ' 

integrating the T dependence of d{e(T))/T. Using the data of figure |l] we obtain the curves 
shown in figure |^ (b) . The line is the result valid for large ./V obtained from the quadratic best 
fit of s c (e). The agreement for N = 300 and 1000 is rather good. We note that increasing TV 
reduces the IS energy range explored by the system for a given fixed Montecarlo simulation 
length, and this in turn reduces the temperature range in figure |^ (b). In general to have 
a good resolution of the IS-energy range N should not be too large, however it cannot be 
too small to have a good sampling of the relevant IS. The sizes reported here give reasonable 
results, even if TV = 48 and N = 1000 are close to the border. 

We finally note that the analysis of the IS done in this Letter can be also useful for systems 
with different types of transition, e.g. the spin glass transition. In figure [| it is shown the 
average IS energy density (e(T)) as function of T for different system sizes for the mean-field 
±J Sherrington-Kirkpatrick model, and in the inset the configurational entropy density s c (e). 
For a more detailed IS-analysis for this model see Ref. |p5| . 

To summarize, in this Letter, using an IS-based analysis, we have shown that finite-size 
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mean-field spin-glass models with discontinuous transition are good candidates for studying 
the glass transition. We note that although the IS approach may be generally applied 
to any systems its relevance for thermodynamics is "a priori" not obvious. For 
example simple one-dimensional models with non trivial configurational entropy 
may show uninteresting thermodynamical behaviour ^7|. Moreover, Stillinger p8{ 
has shown that for realistic systems an IS-based analysis would lead to a zero 
temperature Kauzmann transition. Both problems seem to be absent for the 
finite size mean-field spin glass models discussed in this Letter. 

Our study opens the way for addressing activated processes in glasses through finite-size 
corrections, i.e., calculations beyond the saddle-point approximation, in mean-field spin-glass 
models. Finite-size mean-field spin glass models, for instance p-spin models in the spherical 
approximation, have the double advantage of being analytically tractable for N — > oo and 
easily simulated numerically for finite N, offering a good model to analyze the glass transition. 
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